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Introduction 

The objectives of this paper are two-fold. First, a numerical scheme for the simulation of a 
buoyant, reacting jet is presented with special attention given to boundary conditions. In the 
absence of coflow, a jet flame is particularly sensitive to boundary conditions enforced upon 
the computational domain. However, careful consideration of proper boundary conditions 
can minimize their effect upon the overall simulation. 

Second, results of some preliminary simulations are presented over a range of Froude and 
Damkohler numbers. This range was chosen so as to produce lifted flames in both normal 
gravity and microgravity environments. 

Numerical Method 

A low Mach number approximation [1] is applied to the Navier-Stokes equations and the 
resulting system is solved numerically using a predictor-corrector method similar to that 
described by Najm et al. [2], This predictor-corrector scheme handles large density ratios 
necessary to study buoyancy effects. A one-step, reversible, Arrhenius-type reaction is used 
to model the chemistry. Also, the numerical method is adapted for numerical solution on 
a variable spaced, staggered, cylindrical mesh, with a computational domain as shown in 
Figure 1. For more details about this numerical method, please refer to [3]. 

Boundary Conditions 

Inlet 

The inlet condition corresponds to a jet issuing from a small orifice in a wall. The veloc- 
ity is specified at all points using a “top-hat” profile constructed from a tanh function as 
suggested by Michalke [4], 

Lateral 

As a jet flow develops, it entrains ambient fluid from its surroundings. In the absence of 
coflow, a closed lateral boundary condition encourages recirculation by preventing entrain- 
ment across the lateral boundary. Recirculation of heat and species can have a significant 
effect upon a reacting system. In order to allow entrainment across the lateral boundary, 
a modification of the traction-free boundary condition discussed by Boersma et al. [5] is 
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Figure 1. Diagram of the computational domain with a traction-free lateral boundary con- 
dition and a free-slip collar around the outlet. 


employed. This condition specifies that 


r • n = 0, 


(1) 


where n is the unit outward normal vector at the lateral boundary. The left hand side of 
(1) is proportional to the force due to viscosity exerted on a small surface element of the 
boundary. 


Outlet 

At the outlet of the domain a convective boundary condition was employed of the form: 


dpu 

dt 


+ U„ 


dpu 

dx 


0 


(2) 


The convective speed U max was chosen to be the maximum measured velocity in the plane 
of the outlet [6]. This choice was found to be representative of the small central portion of 
the domain containing the jet with which the simulation is primarily concerned. Therefore, 
it exerts the least amount of unphysical influence upon the region of interest. 


Mass Conservation 

Mass conservation was preserved by correcting pu x at the outlet. The amount of excess 
mass flux j across the outlet was determined by a control volume analysis of the form: 

>=niirt dv+ ih- nds < 3 > 

The outlet mass flux was corrected by uniformly subtracting j /A from pu x , where A is the 
cross sectional area of the outlet. 
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Free-slip Collar 

Physically, the outlet mass flux correction is equivalent to introducing a uniform pressure 
gradient in the x direction across the entire outlet, so that Eqn. (2) can be written: 


dpu 

dt 


+ u ^ 

“ max 

OX 



(4) 


Note that the addition of a corrective pressure gradient is not consistent with the traction 
free boundary condition at the edge of the domain. It was found that the introduction of a 
free-slip “collar” around the outlet (see Figure 1) significantly improves numerical stability 
by decoupling the convective outflow and the lateral traction- free boundary conditions. 

Results and Discussion 

Sixteen axisymmetric simulations were performed with Damkdhler numbers varying in 
increments of 100.0 from 600.0 to 900.0 and inverse Froude number varying in increments of 
0.1 from 0.0 (non-buoyant) to 0.3 (highly buoyant). A mesh of 256 by 512 grid points was 
used with spatial resolution such that in the region of the flame A r ~ 0.03 and Ax ~ 0.04. 
A time step of At = 0.002 was used. Initially, the flame was piloted by the introduction of 
a region of high temperature at the edge of the nozzle. Then, at non-dimensional time 16.0, 
the pilot was turned off, allowing the flame to lift. 

Figure 2 shows contours of density of four resulting lifted flames taken at time t = 192.0 
corresponding to the four different Froude numbers tested and a Damkohler number of 800.0. 
It is immediately apparent that in the case of Fr = 3.33, the flow is much more complicated 
than in the non-buoyant case. Although the flow is unforced at the inlet, the buoyant cases 
develop an instability in the convection layer at the outer edge of the flame. Note that in 
this region, the baroclinic torque has the same sign as the vorticity generation due to viscous 
shear. The initial source of this instability may be numerical noise associated with roundoff 
as well as rather high frequency, small fluctuations in pressure at the outflow boundary. 

Although the non-buoyant momentum-dominated jet is expected to be unstable, it is 
evident from Figure 2 that the instability is not manifested. Although there is always 
numerical noise associated with computations performed on a discrete computer, this noise 
is very low level (on the order of 10 -6 for single precision arithmetic) compared to the 
main flow, so that these instabilities are not triggered. However, the instabilities generated 
by buoyancy seem to be very sensitive to low level excitation. Bahadori et al. [7] note 
that even slightly buoyant flames, with Froude numbers in excess of 40.0, show significant 
differences from non-buoyant flames. This seems reasonable in light of the fact that in a 
laboratory setting low-level noise is usually several orders of magnitude greater than that 
introduced by a computer simulation. Also, the fact that the flow is constrained to be 
axisymmetric in the computer simulations may help stabilize the flow for moderate Froude 
numbers where unstable flow is observed experimentally. Buoyant instabilities are often 
highly three dimensional and this may play a key role in the transition to turbulence for 
lifted flames. 

Conclusions 

From numerical experiments, it was found that a convective outflow boundary condition 
with a choice of U max for the convective wind speed, and a traction-free lateral boundary 
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Figure 2. Density field of lifted flames at time 192 with Damkohler number 800. From left 
to right, the Froude number in each plot is oo, 10.0, 5.0, and 3.33 respectively. 


condition yielded the best performance for the jet flow. The addition of a free-slip collar 
around the outlet provided the decoupling necessary for these conditions to function together. 

It was also observed that the effects of buoyancy tend to destabilize the jet flow. However, 
for moderate Froude numbers, the computer simulated flames appeared to be weakly unstable 
whereas corresponding flames in the laboratory exhibited stronger instabilities. 
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